UMD-PP-94-117 
February 7, 2008 



Predictions for Proton Life-Time in Minimal 
Non-Supersymmetric SO(10) Models: An Updat^\ 



Dae-Gyu Lee and R. N. Mohapatra 

Department of Physics, University of Maryland 
College Park, Maryland 20742 
and 

M. K. Parida and Merostar Rani 

Physics Department, North Eastern Hill University 
P.O.Box 21, Laitumkharh, Shillong-793003, India 



Abstract 

We present our best estimates of the uncertainties due to heavy particle 
threshold corrections on the unification scale M.jj, intermediate scale M/, and 
coupling constant a\j in the minimal non-supersymmetric SO(10) models. 
Using recent CERN e + e~ collider LEP data on sin 2 #vy and a s t r0 ng to obtain 
the two- loop-level predictions for Mu and ajj, we update the predictions for 
proton life-time in minimal non-supersymmetric SO(10) models. 
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1 Introduction 



The hypothesis of a single unified gauge symmetry describing all forces and 
matter at very short distances is a very attractive one from practical as well as 
aesthetic point of view. Right now there are several good reasons to think that 
this gauge symmetry may indeed be SO(10)(l|. The most compelling argument 
in favor of SO(10) comes from ways to explain|| the observed deficit || of the 
solar neutrino flux compared to the predictions @] of the standard solar model in 
terms of a two flavor MSW0 neutrino oscillation. Consistent understanding of the 
data from all four experiments using the MSW oscillation hypothesis requires the 
neutrino masses and mixings to lie in a very narrow range of values. It was shown 
in a recent paper||, that the minimal SO(10) theory that implements the see-saw 
mechanism 0] is a completely predictive theory in the neutrino sector and predicts 
masses and mixings between v e and that are in this range. In addition to this, 
there are many other highly desirable features of the SO(10) theory, such as fermion 
unification into a single {16}-representation, a simple picture of baryogenesis[|[, 
asymptotic parity conservation of all interactions, etc. In view of these, we have 
undertaken a detailed quantitative analysis of the symmetry breaking scales of this 
minimal model in order to pinpoint its predictions for proton life time, especially 
the unertainties in it arising from unknown Higgs masses in the theory. 

Since we are going to discuss the minimal SO(10) model, let us explain what 
we mean by the word 'minimal'. It stands for the fact that a) the Higgs sector is 
chosen to consist of the smallest number of multiplets of SO (10) that is required 
for symmetry breaking and b) only those fine-tunings of the parameters needed to 
achieve the desired gauge hierarchy are imposed. The above fixes the Higgs mass 
spectrum of the model completely. 

Before we proceed further, we wish to make the following important remark 
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about the minimal SO(10) model. For a long time, it was thought that this model 
cannot be realistic since it predicts the relations among fermion masses such as 
m s = and rrid = m e at the GUT scale Mjj, which after extrapolation to the 
weak scale, are in complete disagreement with experiment. However, it was shown 
in Ref. || that in the minimal SO(10) models where the small neutrino masses 
arise from the see-saw mechanism |7]], there are additional contributions to charged 
fermion masses, that solve this problem. They arise from the fact that the (2,2,15) 
submultiplet of the {126}-dim. Higgs multiplet used in implementing the see-saw 
mechanism automatically acquires an induced VEV without additional fine-tuning. 
These additional contributions correct the above mass relations in such a way as 
to restore agreement with observations. The same theory, as mentioned above, also 
predicts the interesting values for neutrino masses and mixings making the minimal 
SO (10) models not only completely realistic but also testable by neutrino oscillation 
experiments to be carried out soon. 

Next, let us mention a word on our choice of non-supersymmetric version of the 
model. While the question of gauge hierarchy certainly prefers a supersymmetric 
(SUSY) SO(10) model, in the absence of any evidence of supersymmetry at low 
energies as well as for the sake of simplicity alone, we believe that minimal non- 
SUSY SO(10) should be thoroughly explored and confronted with experiments. 

Another interesting point that needs to be emphasized is that for the minimal 
set of Higgs multiplets, SO(10) automatically breaks to the standard model via only 
one intermediate stage, that consists of the left-right symmetric gauge group with or 
without the parity symmetry ||10|| , depending on the Higgs multiplet chosen to break 
SO (10). This leads to the following four possibilities for the intermediate gauge 
symmetry: 

A) G 224D = SU(2) L x SU(2) fl x SU(4) c x D 

B) G 224 = SU(2) L x SU(2) R x SU(4) C 
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C) G 2213D = SU(2) L x SU(2) R x U(l) B _ L x SU(3) c x D 

D) G 2213 = SU(2) L x SU(2) fl x U(l) B _ L x SU(3) C 

Case A arises if the Higgs multiplet used to break is a single {54}-dimensional 
onellTJ. Cases B and C arise if a single {210}-Higgs multiplet is used. Depending 
on the range of the parameters in the Higgs potential, either case B or case C arises as 
the intermediate symmetry |12||. Case D arises when one uses a combination of {45}- 



and {54}-dimensional Higgs multiplets|T3|]. The rest of the symmetry breaking 
is implemented by a single {126}-dimensional representation to break SU(2)rX 
XJ(1)b-l as well as to understand neutrino masses and a single complex {10} to 
break the electroweak SXJ(2) L x U(l)y down to U(l) em . These four cases therefore 
represent the four simplest and completely realistic minimal SO(10) models. In the 
rest of the paper, we present calculations of the predictions for proton life-time (r p ) 
in these models as well as the uncertainties in these predictions due to unknown 
Higgs masses and the uncertainties in the low energy input parameters, in order to 
see if the next round of proton decay search at Super-Kamiokande (SKAM)||14|| can 
test this model. 

2 Computation of the Threshold Uncertainties in 
Mu and Mj 

The two main equations in our discussion are i) the two-loop renormalization 
group equation for the evolution of the gauge coupling , i.e., 



do; i &i 2 , 2 r-t \ 



and ii) the matching formula at the mass scale where the low energy symmetry 
group enlarges [IH 
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In Eqs. (1) and (2), oti is the "fine-structure" constant corresponding to the 
gauge group G« and 



Af = Tr^ 2 + Tr^f 2 ln^, (3) 

where Qf is the representation of the gauge group in the representation of the 
Higgs submultiplet H. The expressions for a« and bjj for the four cases are given in 
Table [17j . In deriving the values of a^ and by in various cases as well as to obtain 
the threshold corrections Aj, we need to know the order of magnitude of the mass 
of the various Higgs submultiplets in the models. We obtain these by invoking the 
survival hypothesis for the Higgs multiplets as dictated by the minimal fine tuning 
condition for gauge symmetry breaking ||18||. Using this hypothesis, in Tables Ila-IId, 
we list the various Higgs multiplets whose masses are near the relevant symmetry 
scales along with their contributions to Af. 

We proceed as follows: first using the two-loop equation, we derive the mean 
values for the mass scales in various cases. These results already exist in the 



literature [|T6| , |TT| , I5[ |20| based on the earlier LEP results. In Table III, we have pre- 
sented their values from Ref. which uses the inputs a x {M z ) = 0.016887± 0.000040 
a 2 {M z ) = 0.03322±0.00025; a 3 {M z ) = 0.120±0.007, for further use in calcuating 
r p . These values of Mu and Mj were obtained using analytic integration of Eq. (1) 
which has been done exactly for case A. For cases B, C, D, we have ignored terms 
whose effect in the final result of the renormalization group equation is smaller than 
the error coming from low energy LEP data by a factor of ten or more. We have also 



checked that inputting the most recent LEP |2T| gives results for the mass scales 
which are within the level of accuracy of our calculations. For instance, for Mjj 



the changes are 10' 08 , 10' 01 , 10' 06 , and 10 14 for cases A, B, C and D respectively. 
Then, we estimate the uncertainties in Mj and Mjj due to both the experimental 
uncertainties in the low energy parameters a s , sin 2 #vy and a em as well as the un- 



known Higgs boson masses. For the cases B and D, these were discussed in Ref. |22 ] 
although we refine these uncertainties somewhat, but the threshold uncertainties for 
cases A and C are new. These uncertainties are presented in Table IV. We have 
allowed the Higgs masses to be between 1/10 to 10 times the scale of the relevant 
symmetry breaking . 

The formulas for the threshold effect on the mass scales Mj and Mjj, which 
were used to obtain Table IV, are given below for each symmetry breaking chain. 
We have defined rji = In (Mj^/M), where M is the relevant gauge symmetry breaking 
scale near M^. 
Model A 



A\n(M c+ /M z ) = 077x0 + 0r] 126 + 0r) 54 

+0% + 0.431818% - 0.454545% 
A\n(Mu/M z ) = -0.04r/ 10 - 0.08r/ 126 - 0.04r/ 54 

+0.02^ - 0.106818% + 0.194545% (4A) 

Model B 



Aln(M c /M z ) = +0.0518135% + 0.103627% 6 + 0.0518135% 

-0.025906777^ + 1.12953% - 1.29534% L 
Aln(Mu/M z ) = -0.0621762t7 10 - 0.124352t7 126 - 0.0621762/7210 

+0.0310881% - 0.40544% + 0.554404%,, (4B) 

Model C 
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A\n(M R+ /M z ) = +0.0952381r7io - 0.0952381r7 126 + 0r} 210 

-0.047619^ + 0. 190476^ - 0.142857r7 L1 
A]n(Mu/M z ) = -0.0544218r/io- 0.159864r7 126 - 0.0357143r/ 2 io 

+0.0272109^ + 0.0340136r7 m + 0.117347r7 L1 (4C) 

Model D 



A\n(M R /M z ) = +0.124138r7 10 - 0.0827586r7i2 6 + 0.00689655?745 

-0.062069?^ + 0.22069r7 m - 0.193103% A 
Aln(Mu/M z ) = -0.0781609r7 10 - 0.170115r7 126 - 0.0413793r7 45 

+0.0390805% + 0.0091954%! + 0.158621% A (4D) 

In obtaining the above equations, we have assumed that the particles from a 
single SO(10) representation which have masses of the same order are degenerate. 



This is the same assumption as in Ref . |22fl . Before proceeding to give our predictions 



for proton life-time, few comments are in order. 

a) We want to clarify how we get the uncertainties presented in Table IV. First, as 
already mentioned, we chose M#/Mj or M#/M[/ to vary between 10 _1 and 10 +1 . 
The maximum values of the uncertainties are obtained by allowing the different r/'s 
to vary independently to their extreme values that lead to the largest positive or 
negative uncertainty. The only exception to this is the two parameters r?A L and 
fjH A , which are always kept negative. (See below.) Secondly, in chains A and C 
we do not assume that the left- and right-handed Higgs submultiplets to have same 
mass (that would lead to t]l = i]r)- The reason is that since the masses are close to 
the intermediate scale, where left-right symmetry is broken, the multiplets need not 
necessarily be degenerate. If we assumed the degeneracy, there would be cancellation 
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between r\y and rjji terms reducing the threshold uncertainties [23]; the uncertainties 
we present in table IV are, therefore, most conservative. 

b) In cases B and D, since D-parity is broken at the GUT scale, the masses of At 
in Table Ilb-(l) and Ha in Table Ild-(l) are always above the scale M/, but below 



M[/|22j. Although a priori M& L could be bigger than My, we have kept it smaller 
in presenting the uncertainty in r v . Therefore, r]^ L and t]y A in Eqs.(4B) and (4D) 
are always negative, since we use M = My to define them. In any case, from an 
experimental point of view, the upper value of the uncertainty is not too relevant. 

c) The first set of entries in Table IV is obtained by maximizing the uncertainty in 
Mj whereas the second set is obtained by reversing this procedure. 

d) Note that, in case A, the intermediate mass scale Mj and the unification scale 
My are so close that one might think of this as an almost single step breaking. This 



is similar to the D-parity broken scenario (case B) recently discussed in Ref. p4 
For proton life-time estimate, this is inconsequential. 



3 Predictions for Proton Life-Time 

Now, we present our predictions for proton life-time in the four SO(10) models 
A - D. For this purpose, we need the values of and ay and remember that in 
SO(10) there are extra gauge bosons contributing to proton decay compared to the 
SU(5) model. We use the following formula from the review by Langacker|[25||, where 
the original literature can be found. We write 

r =r (0) F 
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where F p denotes the uncertainty arising from threshold corrections as well as the 
experimental errors in a s , a em , and sin 2 9 w . From Ref. [p5|| , we get for 



Including the F p -factors, we present below the predictions for proton life-time 
in SO(10) (noting that a^j U ^ ~ c%°^ 10 ' ) ). The first uncertainty in the predictions 
below arises from the proton decay matrix element evaluation whereas the second 
and the third ones come from LEP data and threshold correction, respectively [f26| 
Model A 



T p ^ e+n o = 1.44 X l 32-l±-T±1.0±1.9 yrs 



Model B 



r p ^ e+n o = 1.44 x 10 37 - 4± - 7±1 '°- 5 5 oyrs. 



Model C 



Tmp+t0 = 1.44 x io 34 - 2± - 7± - 8±1 - 7 yrs. 



1 p— >e" r 7r' J 



Model D 



Tp ^ e+7r o = 1.44 x 10 37 - 7± - 7± - 9 -2 7 oyrs. 



4 Conclusion 

In conclusion, we have computed the threshold uncertainties in both the inter- 
mediate and the unification scales for all four possible minimal non-supersymmetric 
SO(10) models A - D. We then update the predictions for proton life-time in all 

these cases including the most conservative estimates for the threshold uncertainties 
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in it. We see that for case A, t p is very much within the range of Super-Kamiokande 
search even without threshold corrections. On the other hand, for cases B and C, 
the threshold uncertainties have the effect of bring it within the range of SKAM 
search. 
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Table Caption 

Table I: One- and two-loop /^-function coefficients for models A - D. 

Table Ha: The heavy Higgs content of the model A. The G224 submultiplets in (1) 
acquire masses when SO (10) is broken, while the G123 submultiplets in (2) become 
massive when G224D is broken. The multiplets <f), Rj, and L, in (2) arise from 0(2,2,0) 
in {10}, from A fi (l,3,10), and A i (3,l,10) in {126}, respectively. Also listed, in the 
extreme right column of the tables, are the threshold contributions Aj of the different 
multiplets. 

Table lib: The heavy Higgs particles in the model B whose intermediate symmetry 
is G224- The particles whose masses are on the order of M;y are listed in (1), and 
the particles on the order of Mj are listed in (2). Also listed are their threshold 
contributions. 

Table lie: All the heavy Higgs particles in the model C whose intermediate symmetry 
is G2213D- The submultiplets with masses of order Mjj are presented in (1), and the 
particles with masses of order Mj are listed in (2). The entries in the extreme right 
column denote their threshold contributions. 

Table lid: All the heavy Higgs particles in the model D whose intermediate symme- 
try is G2213- The submultiplets with masses of order Mu are presented in (1), and 
the particles with masses of order Mj are listed in (2). The entries in the far right 
column give their Aj contributions. 

Table III: The values of Mu, Mj and a\j obtained by solving the two- loop renor- 
malization group equation Eq. (1) for models A - D. The results were taken from 



Ref. p0|. 



Table IV: The threshold uncertainties due to the difference between symmetry break- 
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ing scale and masses of Higgs bosons on the order of that scale. For the cases where 
threshold effects in Mj and M v are maximized, corresponding threshold uncertain- 
ties are given in the first two lines and the last two lines, respectively. 
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Table I. 



Model a; 

A ill 11 —Hi 1 765 -j rq 584 7651 rl53 153 1759 1 1 

13)3531 II 3 ' °> 2 J ' «- ' 3 ' 2 J ' t- 2 ' 2 ' 6 J J 

R r 11 23 1 J Jo 45 1 Jo 584 765 \ r9 153 643 \\ 

C {— |, — 1,7, — 7} {{y , 3, y, 12}, {3, y , y , 12}, {y, y, -y , 4}, {|, |, y — 26}} 

D {-3,-|,y,-7} {{8,3,|,12},{3,y,y,12},{|,y,y,4},{|,|,i,-26}} 



Table Ila-(l). 



SO(10) representation G 2 24 submultiplet {A^A^A^} 



10 


H(l,l,6) 


{ 0, 0, 2} 


126 


Co(2,2,15) 


{30,30,32} 




S(l,l,6) 


{ 0, 0, 2} 


54 


S s (3,3,l) 


{ 6, 6, 0} 




S ( (l,l,20') 


{ 0, 0, 8} 




s+(i,i,i) 


{ 0, 0, 0} 
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Table IIa-(2). 



SO(10) representation G123 submultiplet {^[y ^ii^ic} 



1 ( 1 




r 3 1 nl 
{5' °) 


Ri(-2 A 


A3" 1 1 \ 
/5' 1 ' 1 ) 


r 24 n nl 






r 2 n 11 

{5' °' x l 


T-) / 4 


/~3 1 Q\ 


r 32 n 11 

0, 1} 


R 4 (— |a 




{f,0, 5} 


M+h 




{f,0, 5} 




/f.1.6) 


{f,0, 5} 


Li(+V 


f.3,1) 


{f,4, 0} 




/f,3,3) 


{§,12, 3} 




/§A6) 


{f ,24, 15} 



Table Ilb-(l). 



SO(10) representation G 2 24 submultiplet {A^A^A^} 



10 


H(l,l,6) 


{ 0, 0, 2} 


126 


Co(2,2,15) 


{30,30,32} 




S(l,l,6) 


{ 0, 0, 2} 




A L (3,l,l0) 


{40, 0,18} 


210 


S L (3,1,15) 


{30, 0,12} 




S fl (l,3,15) 


{ 0,30,12} 




Ci(2,2,10) 


{10,10,12} 




C 2 (2,2,T0) 


{10,10,12} 




C 3 (l,l,15) 


{ 0, 0, 4} 




s'(i,i,i) 


{ 0, 0, 0} 
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Table IIb-(2). 



G123 submultiplet 
0, Ri, R 2 , R 3 , R 4 , R 5 , Re in TableIIa-(2) 
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Tabic IIc-(l). 



SO(10) representation G2213 submultiplct {^2L^2R^ix^3c} 
10 



126 



210 



Ti(l,l,+iJ|, 3) 


{ 0, 0, 1, 


1} 


T 2 (l,l,-k/|, 3) 


{ 0, 0, 1, 


1} 


Hiij(l,3,+i A /|, 6) 


{ 0,24, 6, 


15} 


Hi /" f 3.1 . — 77 a / 77 .6) 


{24, 0, 6, 


151 


v ' ' 3 y 2 ' / 


{ 0, 4, 3, 


3} 


HoT (3,1, + 77* / 77, 3) 


{ 4, 0, 3, 


3} 


H 3 (2,2,+| A /i, 3) 

v ? 5 1 3 \/ 2' / 


{ 6, 6,16 


4} 

J 


EU2 2 --i/^3) 


{ 6, 6,16 


41 


EU(2 2 8) 


{16,16, 


241 


H fi (2, 2,0,1) 


{ 2, 2, 0, 


0) 


H 7 (l,l,+I,/|, 3) 

f V ' ' 1 3 \ / 2' / 


{ 0, 0, 1, 


1} 


Hfifl.l, — 77 A / 77 « 3) 
O \ 5 " 3 \/ 2' / 


{ 0, 0, 1, 


1} 


B L1 (3, 1,0,8) 


{16, 0, 


91 

J 


Bri (1 3 8) 


{ 0,16, 


91 


B T 9 ( 3, 1 , ~~ 77 A / 77 . 3) 

LZ V; 1 3 y 2' / 


{ 6, 0, 6 


2 J 


Bro(1,3, — 77A/77, 3) 


{ 0, 6, 6 


2 J 


B T "3 (3. 1. + 77 A / 77 . 3) 

' 1 3 y 2' / 


{ 6, 0, 6 


11 
2 j 


Bp^(l,3,+ 7T A / 77, 3) 


{ 0, 6, 6 


2 J 


Br 4 (3, 1,0,1) 


{ 2, 0, 0, 


0} 

J 


Br 4 (1,3,0,1) 


{ 0, 2, 0, 


0} 


B 5 (2,2,+i A /i,6) 

v ' ' 3 \/ 2 ' 


{ 6, 6, 4, 


10} 


B 6 (2,2,--U/f,6) 

u V ' ' 3\/ 2 ' 


{ 6, 6, 4, 


10} 

J 


"D /o O 1 /3 Q\ 

*MAA —3 y 


f Q Q O 

{ 3, 3, 2, 




B 8 (2,2,+i v /f,3) 


{ 3, 3, 2, 


2} 


B 9 (2,2 1) 


{ 1, 1, 6, 


0} 


B 10 (2,2, +v /|, 1) 


{ 1, 1, 6, 


0} 


Bn(l, 1,0,8) 


{ 0, 0, 0, 


3} 


B_(1,1,0,1) 


{ 0, 0, 0, 


0} 


B+(1,1,0,1) 


{ 0, 0, 0, 


0} 
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Table IIc-(2). 



G123 submultiplet 
0, Ri, Li in Table IIa-(2) 



Table Ild-(l). 

SO(10) representation G2213 submultiplet {^L^m^ixi^ic} 
10 Ti, T 2 

126 All H's in Table IIc-1 

Ha(3,1,^/§,1) 
45 Si(l,l,0,8) 
S 2 (3,l,0,l) 
S 3 (l,3,0,l) 
S_(1,1,0,1) 



Table IId-(2). 

G123 submultiplet 
0, Rx in Table IIa-(2) 



Table III. 



Model Mj (GeV) M v (GeV) 



A 


10 13.64 


igl5.02±.25 


40.76 ±.16 


B 


10 10.70 


igl6.26±.25 


46.35 i;l 


C 


10 10.16 


igl5.55±.20 


43.86 ±.18 


D 


10 9.08 


10 16 - 42 ±:S 


46.12 tf 6 



{ 4, 0, 9, 0} 
{ 0, 0, 0, 3} 
{ 2, 0, 0, 0} 
{ 0, 2, 0, 0} 
{ 0, 0, 0, 0} 
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Table IV. 



Threshold Uncertainty 


Model A 


Model B 


Model C 


Model D 


M//M? 


-^Q±0.886 


+2.658 

10-0067 


-j_Q±0.571 


+0.690 
10-0.303 




■^g±0.481 


+0.131 
10-1-240 


■^g±0.031 


-0.138 
10-O.179 


M//M? 


-^q±0.886 


+2.658 
10-0067 


■^giO.000 


+0.303 
10+0.083 




■^q±0.481 


+0.131 
1Q-1-240 


-^q±0.429 


+0.179 
1Q-0.497 
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